Hydrogeophysical modelling of Hisarcik (Kütahya) geothermal field,western Turkey     

Serkan Üner  Gülçin Özürlan Ağaçgözgü  Doğa Düşünür Doğan 《Geophysical Prospecting》2019,67(8):2176-2195

Western Anatolia hosts many low-to-moderate and high-temperature geothermal sources in which active faults play a dominant role to control the recharge and the discharge of geothermal fluid. In this study, we used the two-dimensional geoelectric structure of Kütahya Hisarcık geothermal field, and created a conceptual hydrogeophysical model that includes faults, real topographical variations and geological units. The temperature distribution and fluid flow pattern are also investigated. The depth extension of Hisarcık Fault, electrical basement and low resistivity anomalies related to the presence of geothermal fluid are determined by using resistivity studies in the area. Numerical simulations suggest that Hisarcık fault functioning as a fluid conduit primarily enables hot fluid to be transported from depth to the surface. It is shown that the locations of predicted outflow vents coincide with those of hot springs in the area.  相似文献   

空间波数混合域磁异常场积分解三维数值模拟          下载免费PDF全文

李昆  戴世坤  陈轻蕊  张钱江  赵东东  王顺国  凌嘉宣 《地球物理学报》2019,62(11):4437-4450

本文提出一种空间波数混合域磁异常场三维数值模拟方法.该方法利用磁位三维空间域积分为卷积的特点，沿水平方向进行二维傅里叶变换，把空间域磁位满足的三维积分问题转化为不同波数之间相互独立的垂向一维积分问题.保留垂向为空间域，优势之一在于便于浅层单元剖分可适当加密，随着深度增加，单元剖分适当稀疏，可以准确模拟任意复杂地形和磁性体的磁异常，兼顾了计算精度与计算效率；优势之二在于一维积分垂向可离散为多个单元积分之和，每个单元采用二次形函数表征磁化强度，可得出单元积分的解析表达式，计算精度高、效率高.该方法充分利用一维形函数积分的高效和高精度、快速傅里叶变换的高效性及算法高度并行性，实现了磁异常场高效、高精度的数值模拟.设计棱柱体模型，将模型解析解与空间波数混合域法的数值解对比，结果表明该方法计算精度高、效率高.设计了组合棱柱体复杂模型，对比分析了标准FFT扩边法与Gauss-FFT法的计算精度与计算效率，总结了标准FFT的扩边系数选取策略.针对任意复杂地形条件下的磁异常模拟问题，本文提出一种适用于起伏地形条件下的磁异常场快速计算方法，并对其有效性进行了验证.  相似文献   

Travel time computations using a compact eikonal equation for vertical transverse isotropic media          下载免费PDF全文

Peyman Moghaddam  Reza Khajavi  Henk Keers 《Geophysical Prospecting》2018,66(8):1475-1484

Eikonal solvers often have stability problems if the velocity model is mildly heterogeneous. We derive a stable and compact form of the eikonal equation for P‐wave propagation in vertical transverse isotropic media. The obtained formulation is more compact than other formulations and therefore computationally attractive. We implemented ray shooting for this new equation through a Hamiltonian formalism. Ray tracing based on this new equation is tested on both simple as well as more realistic mildly heterogeneous velocity models. We show through examples that the new equation gives travel times that coincide with the travel time picks from wave equation modelling for anisotropic wave propagation.  相似文献   

Response and capacity of monopod caisson foundation under eccentric lateral loads     

Tanmoy Kumar Deb  Baleshwar Singh 《Marine Georesources & Geotechnology》2018,36(4):452-464

Monopod caisson foundation is a viable alternative for supporting offshore wind turbines located at shallow water depths. This foundation system has to resist overturning moment generated due to resultant lateral load, arising from wind and water wave action, that can act at any loading height above the seabed. This paper presents results of a numerical investigation performed to determine the influence of loading height, caisson geometry and superstructure load on the ultimate lateral capacity, initial stiffness, and soil failure zone of the foundation, when installed in very dense sand. Both the ultimate and serviceable states of the caisson foundation obtained from the analyses are represented in terms of envelopes plotted between lateral load and overturning moment. Simplified expressions, which take into account the influence of caisson geometry, loading height, and soil properties, are also presented to serve as a preliminary base for design of the monopod caisson foundation.  相似文献   

Numerical study for vegetation effects on coastal wave propagation by using nonlinear Boussinesq model     

《Applied Ocean Research》2018

The vegetation has important impacts on coastal wave propagation. In the paper, the sensitivities of coastal wave attenuation due to vegetation to incident wave height, wave period and water depth, as well as vegetation configurations are numerically studied by using the fully nonlinear Boussinesq model. The model is based on the implementation of drag resistances due to vegetation in the fully nonlinear Boussinesq equation where the drag resistance is provided by the Morison’s formulation for rigid structure induced drag stresses. The model is firstly validated by comparing with the experimental results for wave propagation in vegetation zones. Subsequently, the model is used to simulate waves with different height, period propagating on vegetation zones with different water depth and vegetation configurations. The sensitivities of wave attenuation to incident wave height, wave period, water depth, as well as vegetation configurations are investigated based on the numerical results. The numerical results indicate that wave height attenuation due to vegetation is sensitive to incident wave height, wave period, water depth, as well as vegetation configurations, and attenuation ratio of wave height is increased monotonically with increases of incident wave height and decreases of water depth, while it is complex for wave period. Moreover, more vegetation segments can strengthen the interaction of vegetation and wave in a certain range.  相似文献   

Upper Ocean Four-Dimensional Variational Data Assimilation in the Arabian Sea and Bay of Bengal     

Suneet Dwivedi  Atul Srivastava  Alok Kumar Mishra 《Marine Geodesy》2018,41(3):230-257

A regional ocean circulation model with four-dimensional variational data assimilation scheme is configured to study the ocean state of the Indian Ocean region (65°E–95°E; 5°N–20°N) covering the Arabian Sea (AS) and Bay of Bengal (BoB). The state estimation setup uses 10 km horizontal resolution and 5 m vertical resolution in the upper ocean. The in-situ temperature and salinity, satellite-derived observations of sea surface height, and blended (in-situ and satellite-derived) observations of sea surface temperature alongwith their associated uncertainties are used for data assimilation with the regionally configured ocean model. The ocean state estimation is carried out for 61 days (1 June to 31 July 2013). The assimilated fields are closer to observations compared to other global state estimates. The mixed layer depth (MLD) of the region shows deepening during the period of assimilation with AS showing higher MLD compared to the BoB. An empirical forecast equation is derived for the prediction of MLD using the air–sea forcing variables as predictors. The surface and sub-surface (50 m) heat and salt budget tendencies of the region are also investigated. It is found that at the sub-surface, only the advection and diffusion temperature and salt tendencies are important.  相似文献   

Effects of the geometry of two‐dimensional fractures on their hydraulic aperture and on the validity of the local cubic law          下载免费PDF全文

Ahmad Hajjar  Luc Scholtès  Constantin Oltéan  Michel Antoine Buès 《水文研究》2018,32(16):2510-2525

Flow through rough fractures is investigated numerically in order to assess the validity of the local cubic law for different fracture geometries. Two‐dimensional channels with sinusoidal walls having different geometrical properties defined by the aperture, the amplitude, and the wavelength of the walls' corrugations, the corrugations asymmetry, and the phase shift between the two walls are considered to represent different fracture geometries. First, it is analytically shown that the hydraulic aperture clearly deviates from the mean aperture when the walls' roughness, the phase shift, and/or the asymmetry between the fracture walls are relatively high. The continuity and the Navier–Stokes equations are then solved by means of the finite element method and the numerical solutions compared to the theoretical predictions of the local cubic law. Reynolds numbers ranging from 0.066 to 66.66 are investigated so as to focus more particularly on the effect of flow inertial effects on the validity of the local cubic law. For low Reynolds number, typically less than 15, the local cubic law properly describes the fracture flow, especially when the fracture walls have small corrugation amplitudes. For Reynolds numbers higher than 15, the local cubic law is valid under the conditions that the fracture presents a low aspect ratio, small corrugation amplitudes, and a moderate phase lag between its walls.  相似文献   

考虑地形起伏和障碍体破裂的汶川地震强地面运动数值模拟   总被引：2，自引：1，他引：1       下载免费PDF全文

赵由佳  张国宏  单新建  尹昊  屈春燕 《地球物理学报》2018,61(5):1853-1862

2008年5月12日四川汶川地区发生M_W7.9地震，震中位置103.4°E，31.06°N.这次地震造成了以汶川、映秀为中心及其周边地域建筑物的严重破坏和人员的重大伤亡，且因为高山等地形复杂区域抢险救灾的艰巨性，为及时救援造成很大干扰.为更好理解地形因素对于强地面数值模拟结果的影响，建立了包含地形起伏影响及去除地形影响的两类模型.同时，依据震源破裂过程运动学反演结果，建立了包含障碍体破裂过程的震源滑动模型，实现断层分段、空间倾角以及滑移角的动态设定.基于动力学的地震动模拟方法，通过对地震波传播过程的数值计算和后处理分析，模拟由地震激发的区域强地面运动过程.结果显示：（1）强震动台站的断层距对地形效应具有放大或抑制作用，距离断层破裂带越近，地形效应越明显，反之，距离越远，则地形效应越微弱；（2）因为地形高差与障碍体的影响，地震造成的峰值可能出现在震中区域之外；（3）考虑地形影响模型的地表峰值速度（PGV）区域位于汶川与北川附近；而未考虑地形影响模型的PGV区域位于灌县—江油断层的后半段，处清平、安县附近；对汶川地震近实时强地面运动波场的模拟、峰值图谱的圈定及未来大地震强地面运动特征的预测都有重要指示意义.  相似文献   

A one‐dimensional process‐based approach to study reservoir sediment dynamics during management operations          下载免费PDF全文

Lucie Guertault  Benoit Camenen  André Paquier  Christophe Peteuil 《地球表面变化过程与地形》2018,43(2):373-386

A comprehensive understanding of the dynamics of erosion and sedimentation in reservoirs under different management conditions is required to anticipate sedimentation issues and implement effective sediment management strategies. This paper describes a unique approach combining fluvial geomorphology tools and morphodynamic modeling for analyzing the sediment dynamics of an elongated hydropower reservoir subjected to management operations: the Génissiat Reservoir on the Rhône River. Functional sub‐reaches representative of the reservoir morphodynamics were delineated by adapting natural river segmentation methods to elongated reservoirs. The segmentation revealed the link between the spatial and temporal reservoir changes and the variability of longitudinal flow conditions during reservoir management operations. An innovative modeling strategy, incorporating the reservoir segmentation into two sediment transport codes, was implemented to simulate the dynamics of erosion and sedimentation at the reach scale during historic events. One code used a bedload approach, based on the Exner equation with a transport capacity formula, and the other used a suspended load approach based on the advection–dispersion equation. This strategy provided a fair quantification of the dynamics of erosion and sedimentation at the reach scale during different management operations. This study showed that the reservoir morphodynamics is controlled by bedload transport in upper reaches, graded suspended load transport of sand in middle reaches and suspended load transport of fine sediments in lower reaches. Eventually, it allowed a better understanding of the impact of dam management on sediment dynamics. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

基于交错网格有限差分法的断面波叠前时间偏移成像影响因素正演分析          下载免费PDF全文

徐嘉亮  周东红  王玉英  边立恩  吕振宇 《地球物理学报》2018,61(2):733-741

断层和断裂带的有效识别是地震资料解释中的重要环节，断层在地震信号响应中以断面波的形式体现，因此断面波成像的质量关系到断层的精细识别与刻画.本文利用精度较高的交错网格有限差分正演模拟方法对断面波成像的影响因素进行了正演研究，主要正演分析的参数包括采集因素中的电缆长度和采集方向，地质因素中的断层倾角、断距、反射系数，以及处理因素中的偏移方法等几个方面.通过正演论证得出：采用合理的采集参数能够提高断面波的照明度；有效结合地质因素能够提高断面波的解释精度；利用合理的偏移方法能够使断层归位更加准确，断面波有效成像.基于以上结论，对于断面波的精确识别与刻画，应综合采集因素，处理因素及地质因素，只有这样才能提高断层的解释精度，有效减小解释误差.  相似文献